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ABSTRACT 

Context. The profiles of emission lines formed in the corona contain information on the dynamics and the heating of the hot plasma. 
Only recently has data with sufficiently high spectral resolution become available for investigating the details of the profiles of 
emission lines formed well above 10 6 K. These show enhanced emission in the line wings, which has not been understood yet. 
Aims. We study the underlying processes leading to asymmetric line profiles, in particular the responsible plasma flows and line 
broadening mechanisms in a highly filamentary and dynamic atmosphere. 

Methods. Line profiles of Fexv formed at 2.5 MK acquired by the Extreme ultraviolet Imaging Spectrometer (EIS) onboard the 
Hinode solar space observatory are studied using multi Gaussian fits, with emphasis on the resulting line widths and Doppler shifts. 
Results. In the major part of the active region, the spectra are best fit by a narrow line core and a broad minor component. The latter 
contributes some 10% to 20% to the total emission, is about a factor of 2 broader than the core, and shows strong blueshifts of up to 
50km/s, especially in the footpoint regions of the loops. On average, the line width increases from the footpoints to the loop top for 
both components. A component with high upflow speeds can be found also in small restricted areas. 

Conclusions. The coronal structures consist of at least two classes that are not resolved spatially but only spectroscopically and that 
are associated with the line core and the minor component. Because of their huge line width and strong upflows, it is proposed that 
the major part of the heating and the mass supply to the corona is actually located in source regions of the minor component. It might 
be that these are identical to type II spicules. The siphon flows and draining loops seen in the line core component are consistent with 
structures found in a three-dimensional magneto-hydrodynamic (3D MHD) coronal model. Despite the quite different appearance of 
the large active region corona and small network elements seen in transition region lines, both show similar line profile characteristics. 
This indicates that the same processes govern the heating and dynamics of the transition region and the corona. 
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1. Introduction 

Understanding the dynamics of the plasma in the corona of the 
Sun is pivotal for unveiling the processes that govern the heating 
of the outer atmosphere of the Sun to temperatures of millions of 
degrees. Imaging instruments can give vital information on the 
evolving structures. However, they only provide information on 
\ the apparent motions, which do not necessarily have to be real 
plasma motions, and they are broad band in wavelength, thereby 
mixing several spectral lines formed at different temperatures. 
Spectroscopy provides details on the emission line profiles and 
thus supplies crucial information for investigating the thermal 
and dynamic structure of the coronal plasma. 

Various scenarios and models make detailed predictions 
on the spectral profiles to be observed. Only selected aspects 
should be mentioned he r e. The nanoflare-heating models of 
iPatsourakos & Klimchukl (I2006I) based on ID hydrodynamic 
loop models predict high-temperature upflows of hot plasma fol- 
lowing transient heating events (introduced adhoc in the model). 
These would show up as a weak additional component in the 
wing of emission lines formed well above 10 6 K with a blueshift 
of so me lOOkm/s (see their Fig. 3). Based on their observa- 
tions, iDe Pontieu et al.l (I2009I) propose a scenario for the mass 
cycle of the active region corona where hot plasma related to 
type II spicules is propelled upwards to supply mass to the 
corona (which subsequently rains down after cooling). This sce- 
nario would also imply a strongly blueshift e d com ponent in the 
corona l emission l ines. S imilarly, iTu et al.l (I2005I) or, more re- 
cently, iTian et all (l2010h have found evidence of the plasma 



in open field regions near the limb above temperatures of sev- 
eral 100.000 K being accelerated to form the open wind, with 
blueshifts of some 30 km/s seen in lines formed at 2 MK. If (lo- 
cally) open regions on the disk, e.g. adjacent to an active region, 
were to produce a similar outflow, this could show up in the spec- 
tra as excess emission in the blue wing of the emission line, too. 

Line profiles originating in the transition region from the 
chromosphere to the corona have been intensively investigated 
using vacuum ultraviolet spectrometers such as HRTS (High 
Reso lution Telescope and Spectrograph; iBrueckner & Bartoel 
119831) or S UMER (Solar Ultrav iolet Measurements of Emitted 
Radiation; IWilhelm et all Il995l) . These provided a high spec- 
tral resolution of about 20 000 or more allowing them not only 
to derive the shifts of the line centroid and width, but also to 
detect peculiarities in the line profiles. The most noticeable fea- 
tures in the tra nsition region line profiles are found during ex- 
plosive events (Der e et all [l989): two distinct satellites to the 
blue and red indicating a bi-directional flow, most likely induced 
by a small reconnection event. The spatial and temporal evolu- 
tion of these were first described by Innes et al. (1997), confirm- 
ing t he nature of a bi-directional flow. Kjelds eth Moe & Nicolasl 
(1 19771) were the first to show that transiti on region sp ectra 
show enhanced emission in the line wings. iPeted (2000) pre- 



sented evidence that these enhanced wings are predominantly 
found above the chromospheric network and are best fitted by 
a double Gaussian with a narrow line core and a broad second 
(or wing) component. A double Gaussian fit with one broad 
component is even more significant than a triple Gaussian fit 
with two narrow components accounting for the wing emission 
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Fig. 1. Intensity maps in selected emission lines and photospheric magnetogram. The line formation temperatures (peak of ion 
fraction, logT) is given in each panel. The images are plotted on an inverse linear scale, scaled with respect to the median intensity 
in each map, (/), with white representing a value l/R and black R where R is given in the plots in units of (J). The maps shown 
here are recorded co- spatially and co-temporally. The lower right panel shows the radiance in the Fe XV line that is analyzed in this 
paper in depth. The lower left panel shows the MDI magnetogram of the photosphere taken at 19: 15 UT during the EIS scan (linearly 
scaled from — 1.5kG to +1.5 kG. The contour lines show two areas with high- velocity components; i.e., the minor component of 
a free double Gaussian fit shows a blueshift higher than 55km/s when smoothed over 11x11 pixels (cf. Sects. 13. 3 1 and 14.41) . The 
rectangles show areas with average spectra plotted in Fig. and the crosses indic ate the posi t ions of the individual spectra (at 
one spatial pixel) shown in Fig.ffl The dotted rectangle indicates the region used by lHara et al.l (120081) for their average spectrum 
(^85" x 1 10"). The intensities are plotted on a logarithmic scale in Fig. IB.ll to show the patterns in the regions with low emission. 



(iPeter & Brkovicll2003l) . This questions the previous interpreta- 
tion of the win g excess as cau s ed by small counterparts of ex- 
plosive events (iDere & Masonl Il993l) . The width of the broad 
wing component see ms to increa se monotonically with line for- 
mation temperature (iPeterl 120011) . which is consistent with line 
broadening due to magneto-acoustic waves. 

Our knowledge of the detailed spectral profiles of lines 
formed in the hotter corona is more limited by instrumenta- 
tion in the era o f SoH O (Solar and Heliospheric Observatory; 
iDomingo et all 1 19951) and before . The Coronal Diagnostic 
Spectrometer (CDS: lHarrison et all 1 19951) did not provide suf- 
ficient spectral resolution to reveal details in the line pro- 
files (except for extreme cases, e.g. flares), and the SUMER 
spectra of lines formed at hot temperatures are problem- 
atic in terms of signal-to-noise ratio. This situation changed 
with the the Ex tr eme u lt raviolet Imagi n g Spec trometer (EIS; 
iKorendvke etall 120061: ICulhane et all l2007h onboard the 



Hinode solar space observatory (iKosugi et all 12007). While the 
spectral resolution of about 4000 is worse than SUMER or 
HRTS, the good signal-to noise ratio provides the possibility of 
studying the details of the line profile. 

The first study of the asym metries of coro nal lines in an ac- 
tive region was conducted by lHara eta l. (2008) using EIS data. 
As a first approach they performed single Gaussian fits to two 
emission lines and found deviations from a single Gaussian, 
showing excess emission in the blue wing near loop foot 
points. Consequently they int erpret their findings by favorin g 
the nanoflare heating model by lPatsourakos & Klimchukl(l2006l) . 
However, as demonstrated in this paper, a closer inspection of 
the data by a more advanced data processing shows a more de- 
tailed picture. Indeed, evidence is found for upflows near the 
loop footpoints, but these are probably not consistent with the 
nanoflare scenario. 
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Fig. 2. Spatial maps of line shift and width for Fe XV (284 A). These are shown for single Gaussian fits (left column), as well as for 
the core component and the minor (wing) component for a free double Gaussian fit. The panels are co-spatial and co-temporal with 
respect to Fig.Q] The respective color scales for the single Gaussian and the core of the double Gaussian are the same, while the 
scales for the minor component of the double Gaussians span far wider ranges. The Doppler shifts have been corrected for the tilt 
of the spectrograph slit and the orbital motion of the spacecraft, while the line widths (exponential widths) have not b een corr ected 
for instrumental broadening. The crosses, rectangles, and contours have the same meaning as in Fig.Q] See Sects. l3.2l and l3.3l 



The paper is organized as follows. First the data and their 
reduction are discussed in Sect.[2] before the multi Gaussian line 
profile fits are analyzed in Sect.[3] Subsequently, the results are 
discussed with respect to the high-velocity outflows in the ac- 
tive region periphery (Se ct.14.41) , to the heating and mass supply 
in the active region (Sect. l4.5l L and to the self- similarity of pro- 
cesses in the active region corona and the quiet Sun network 
transition region (Sect. 14.61) . Finally, the results are put into the 
context of recent three-dimensional magneto-hydrodynamic (3D 
MHD) coronal models in Sect. 14.71 and the conclusions summa- 
rized in Sect.0 



2. Observations and data processing 



To ensure a good comparison with existing work the present 
study deals with the same data set as analyzed previously by 
Hara_etaL|(|2008), but now with a new more complex method to 
investigate the spectral profiles. A raster scan is analyzed cov- 
ering about 234" x 256" near disk center (cf. Fig.Q]) including 
active region NOAA 10938. The \" wide slit was used with a 
sampling of 1 "/pixel along the slit and a raster step of about \" 
to build the map shown in Figs.Q]and[2l The exposure time was 
30 s. The data were taken on 18 Jan. 2007 starting at 18:12 UT, 
and it took about 135 minutes to complete the raster scan. The 
maps shown in Fig.Q] were shifted along the slit (^-direction) to 
account for the spatial misalignment between the two EIS de- 
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tectors (hence the white stripes at the top of some maps; Kamio 
2O1C0). 

Because the goal of this study is to investigate the line pro- 
files originating from the hot part of the corona, this study will 
concentrate on the Fe XV line at 284 A, formed at about 2.5 MK 
under ionization equilibrium conditions. Some sa mple spectra 
are shown in Figs.[3] and HI As discussed in Sect. 14.11 and ap- 
pendix [A] there are potential blends in the blue wing of the Fe XV 
line by Fe XVII and Al IX. However, under active region condi- 
tions, as well as for quiet Sun or flare conditions, these should 
contribute substantially less than the line wing excess found for 
Fe XV, and these blends are at higher blueshifts than the addi- 
tional spectral component needed to fit the Fe XV line profile. 
Thus the line asymmetries discussed in this paper cannot be at- 
tributed to line blends. 

The raw data were processed using the procedures and cali- 
bration data supplied in SolarSofQ by the instrument team. This 
included removing the CCD pedestal and dark current, flagging 
problematic pixels, as well as calib rating and estimating errors. 
For further details see I Young et al.l (120091) . 

For investigating the line profile, single and double Gaussian 
fits with a continuum have been applied. As the goal is to de- 
tect additional spectral components accounting for small excess 
in the line wings, a reliable fitting procedure has to be used, 
as normal fitting algorithms (e.g. based on hill climbing) may 
fail. Thus t he genetic alg o rithm -based optimization method de- 
scribed by CharbonneauJ (1 19951) is used. This method also al- 
lows constraining the line fit parameters, e.g., to force the line fit 
to find solutions only with one spectral component at high blue 
shifts. The genetic algorithm is a global method, which is why 
it is effective in finding the global minimum of a given function 
(here the difference between line fit and data). This contrasts 
with local methods, which are prone to getting stuck in local 
minima. This is of interest especially if the data to be fitted are 
noisy or the fitti ng problem is not determined very well (see e.g. 
the analysis by Peter, 1999, 2000). The latter is the case here, 
especially when fitting a double Gaussian (seven free parame- 
ters) to a line profile with 24 wavelength pixels across the pro- 
file (cf. Figs. [3] and |4]). To estimate the errors, a local optimiza- 
tion method is used (gradient expansion algorithm as provided 
by the curvefit procedure of IDL, Interactive Data Language). 
Using the results from the genetic algorithm as an initial guess, 
the local method can provide a good fit and automatically returns 
the standart deviation for the fit parameters. 

The Doppler shifts have to be corrected for the line shifts 
induced by the orbital motion of the spacecraft, as well as for 
the inclination of the slit. For this we used the Doppler shift of 
the line core, i.e., of the core component as determined from the 
double-Gaussian fits. First the slit inclination is corrected for ac- 
cording to the SolarSoft routines (by P. Young), then the orbital 
variation is determined by averaging the raster map along the slit 
(^/-direction), by this determining the sinusoidal-like variation in 
time that corresponds to the variation along the ^-direction in the 
raster map. This average line shift along the x-direction (viz. in 
time) is then subtracted from the originally determined Doppler 
maps. In the absence of a possibility for an absolute wavelength 
calibration, the average shift (of the line core) is set to zero. Thus 
the Doppler shifts discussed here are only relative to the average 
field of view. Previous studies of the quiet Sun and active regions 
show that the absolute line shifts at about 2 MK can be expected 
to range from 5 km/s to 10 km/s towards the blue (iPeter & Judgel 

1 Private communication. A publication is in preparation. 

2 http://www.lmsal.com/solarsoft. 



1 19991: iTeriaca et all 1 1999h . Thus the Doppler shifts derived in 
this paper might have to be shifted by this amount to the blue 
to obtain an absolute scaling. However, most of the line shifts 
discussed here will be on the order of 50 km/s, and thus this cor- 
rection might not be of vital importance for the present study. 

3. Multi-Gaussian line profile fits 

When performing a multi-Gaussian fit to find an optimal rep- 
resentation of the observed line profile, one makes an implicit 
assumption, namely that a finite number of spatial regions with 
different plasma properties contribute to the spectrum seen in 
the spatial resolution element. In the case of a double Gaussian, 
it could be that this is "background" providing the main emis- 
sion (the line core), and a small part of the resolution element 
could be filled with a strand of strong upflow (showing up in 
the line wing). In principle, this could also come from temporal 
variation, but this seems unlikely within the 30 s exposure time 
data discussed here. 

Other methods of investigating the line profiles share this 
character. If one does a single G aussian fit and investigates 
the residual (e.g. lHara et all 120081) . one is assuming two spa- 
tial components, too. Also when comparing the emission in 
the red and blue wings to char acterize the line asymmetry, e.g. 
iMcIntosh & De Pontieul d2009) assume that this asymmetry is 
caused by a second spatial component in the resolution element. 

In contrast to this assumption that the emission is composed 
of distinct spatial regions alternative interpretations are also pos- 
sible. The line asymmetries could be produced by opacity ef- 
fects — this is why photospheric lines are often investigated 
using the bi- sector technique. However, the extreme ultraviolet 
lines are optically thin, in general. Only very strong lines un- 
der special conditions show o pacity effects, e .g., the strong C IV 
1548 A line near the limb (e.g. lMari ska. 1992). Another possibil- 
ity would be that the shape of the line profile reflects the veloc- 
ity distribution function of the ions, which is briefly discussed in 
Sect. 14.31 Finally an asymmetric instrumental profile could cause 
the asymmetric line profiles, but as discussed in Sect. 14.21 this is 
not probable for the data under investigation here. 

In the following the properties of a single and a free dou- 
ble Gaussian fit to the data will be investigated first (Sect. 13.21 
and !3.3K a constrained double gaussian fit with one co mpon ent 
forced to be at high blueshifts will be studied later (Sect. 13.41) . 

3.1. A short note on the line widths 

It is assumed that each component of the spectral profile can be 
described by a Gaussian, 

^Gauss (v) = I exp f — VD ^ ] , (1) 

where the wavelength is given in Doppler shift units, v, and the 
free parameters are the peak intensity, position, and width, /, v®, 
and Wi/ e . Please note that this Gaussian width w\/ G is related to 
the full width at half maximum, Fwhm, of the line through 

Fwhm = 2V / In2 wi/ e . (2) 

Following iDoschek et al.l (|2007), the instrumental broaden- 
ing of the EIS spectrometer is about 56 mA Fwhm, which trans- 
lates into a Gaussian width of about 35 km/s for the line under 
investigation, i.e Fe XV at 284 A (cf. Table[T]). The data shown in 
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Fig. 3. Average spectral profiles of Fe XV (284 A). These are averages of the regions shown by small squares in Fig. Q] labeled (a) - 
(d), each covering 10" x 10"(100 individual spectra). The spectra are shown on a logarithmic scale to emphasize the excess in the 
line wings. The top and middle rows show the spectra (diamonds) along with the errors (bars) normalized to the continuum intensity. 
(The spectra in each column are identical, e.g. panels al and a2.) The respective panels in the top row show the double Gaussian 
fit (solid) along with the two components (dashed), the middle column shows the single Gaussian fit (solid). The numbers with the 
plot give the contribution of the respective component to the total line intensity /, the Doppler shift and the exponential width 
Wi/ e . The lower row shows the residual of the respective single Gaussian fits. For regions (a) and (b) these residuals show a roughly 
Gaussian profiles with peaks near — 120 km/s as indicated by the solid lines (panels a3, b3; numbers give the fraction of the residual 
when compared to the total line intensity). See Sects. l3.2l and l3.3l 



the figures in this study are not corrected for instrumental broad- 
ening, because the main aim of the paper is not a quantitative 
interpretation of the derived line widths. In order to characterize 
the line fits and to compare them with the actually observed spec- 
tral profiles, it is preferable to keep the line widths as returned 
from the fits. Where necessary the respective values corrected 
for instrumental broadening are quoted. In principle (assuming 
Gaussian instrumental profiles), the "solar" line width is given 
by (w\ t — w? nstr ) 1//2 , where is the line width returned by the 
fitting procedure and u>i nstr the instrumental width. 

For comparison, TableQ] also lists the thermal line width of 
Fe XV at 284 A, assuming that the line is formed near the temper- 
ature of maximum ion fraction, along with the total line width to 
be expected if thermal broadening would be the only broadening 
mechanism. 



Table 1. Line widths related to Fe XV (284 A). 





FWHM 


exponential width 




[km/s] 


[km/s] 


instrumental width" w- m ^ x 


59 


35 


thermal width 6 ^thermal 


45 


27 


total width (wi sti + w t 2 her mai) 1/2 


75 


45 



a According to lDoschek et all (I2007T) : FWHM = 56 mA; at 284 A. 
b ^thermal = 2ksT / /jl with line formation temperature logT[K]=6.4 
and molecular weight \i = 55.8 amu. 
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Fig. 4. Individual spectral profiles of Fe XV (284 A). These were recorded in one spatial pixel x \") shown by crosses in Fig.Q] 
labeled (a) - (S). As in Fig. [3] the profiles are shown on a logarithmic scale to emphasize the line wings. The top row shows the line 
profiles along with the free double Gaussian fits (cf. Sect. l3.3lh the bottom ro w for the constrained double Gaussian fits where the 
minor component was forced to have a blueshift of at least 70 km/s. (cf. Sect. 13.41) . The dashed lines show the two components of 
the fit, the solid line the sum of the two components. The diamonds show the actual spectra with the bars representing the errors. 
The observed spectra in each column are identical (e.g. the spectr a in al and a2). As in Fig.[3j the numbers in the plots give the line 
fit parameters, now accompanied by the values for Xr- See Sects. 13. 3l and l3.4l 



3.2. Single Gaussian fits and excess emission in the line 
wings 

As a first step, just a single Gaussian fit is preformed. The left- 
most column of Fig.[2] shows the resulting maps in line shift and 
width The maps are very similar to tho se shown in [Hara et al.l 
( 2008), of course. As already noted by lHara et al.1 (120081) . the 
spectral profiles are in general not well fit by a single Gaussian 
and show excess emiss ion in the li ne wings. In the sample (aver- 
age) spectra shown by lHara et al.l (120081) . the line excess seems 
to be concentrated in the blue wing of the line, and they inter- 
pret this as the result of chromospheric evaporation as expected, 
e.g., following impulsive he ating events. This would be expe cted 
according to the models of iPatsourakos & Klimchukl ([2006). 

Th e region used for the average spectrum in the Har a et all 
(2001 ) study is outlined in Figs.Q] and [2] (dotted rectangles). 
While Hara et aD (2008) note that this area represents the foot- 
point region of the active region loop system, it is clear from 
inspection of Fig.Q] that the situation is far more complex. 
Therefore in the present paper smaller areas covering about 
10" x 10" will be examined. They contain 100 spatial pixels and 
are outlined in Fig s.[T]and[2]( small rectangles). Only region (b) is 
covering a loop footpoint. Area (a) is in a fan-like structure seen 
in the cooler lines (Fe VIII and Fe x), while region (c) is located 
in the middle of a hot loop (seen in Fe XV). Area (d) is located 



in the middle of the active region, where the coronal emission 
is only very weak. The important observation is that these small 
areas show quite different line profile characteristics. 

In Fig. [3] average spectra in these selected areas (a) - (d) are 
displayed. The middle row shows the raw spectra along with the 
single Gaussian fits an d the lowe s t row displays the residuals of 
the fits. In contrast to lHara etaD (120081) . the spectra are plotted 
on a logarithmic scale to give a better representation of the line 
wing excess. It is clear that only a fraction of all spectra show 
an excess emission in the blue wing alone. Only the case (a) 
labeled "upflow region" in the fan-like structure shows a clear 
blue- wing-only excess (indicated by the Gaussian for the resid- 
ual; Fig. [5J panel a3). In general, excess emission is found in both 
wings of the line, sometimes with more weight to the blue (as 
in panel b3), mostly with relatively comparable excess in both 
wings (panels c3 and d3), but almost never with an excess in 
only the red wing. 

The asymmetries of the line profile will also influence the 
line shifts and widths derived from the single Gaussian fits, of 
course. For example, there will be the tendency for a line profile 
with excess emission in the line wing to result in a larger line 
width for a single Gaussian fit: the fitting procedure will push 
the line width to values that are as high as possible to account 
for the excess line wing emission, even though it cannot fully 
account for the excess. 
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Fig. 5. Distribution of fit quality and contribu- 
tion to line intensity. The left panel shows the 
histograms for the distribution of the quality 
of the fits Xr to the Fexv line profile. The 
thin solid line shows the theoretical distribution 
(with Xr scaled by 0.5 to account for the prob- 
able error in the EIS error determination). 
The right panel shows the contribution of the 
second minor components of the constrained an 
the free double G aussi an fits to the total line in- 
tensity. See Sects.l331andl431 



The single Gaussian fits show a systematic change in the line 
width from low values near the middle (top) of the hot loop sys- 
tem (near rectangle c in Fig.Q] Fe XV panel, and in Fig. [2 lower 
left panel) to the footpoints of t he loops (re s pectiv e rectangles b). 
This has been already noted by Har a et al.l (120081) . who interpret 
this as evidence o f increased heating near th e loop footpoints as 
advocated e.g. by lAschwanden et all (120071) . 

However, this observational finding of lHara eta! ] (120081) is 
only an artifact of the excess emission in the line wings. When 
performing a double Gaussian fit, the line core component shows 
the opposite trend in the line width, with smaller widths near the 
loop footpoints (cf. lower middle panel of Fig. [2]) ! This does not 
rule out the possibility of the heating being concentrated near 
the footpoints, as some mechanism has to account for the excess 
line wing emission, but it clearly shows that more care has to be 
taken when analyzing EUV emission line spectra with a spectral 
resolution and a signal-to-noise ratio good enough to see line 
asymmetries. 



be comparable. Mostly this is correct, but not in regions where 
the minor component shows strong blueshifts. Here the single 
Gaussian fit shows strong blueshifts, while the double compo- 
nent fit reveals that the line core is mostly redshifted (!), and the 
minor component is a satellite in the blue wing (for sample spec- 
tra, see panels al and al in Figs. [3] and |4]). This strong outflow 
component is restricted to two small areas and is discussed in 
more detail in Sect. l4.4Fl Again, this emphasizes the need for a 
proper multi-Gaussian line fit. 

The line widths of both the core and the minor component 
are larger in the middle part of the active region than in the foot- 
point regions; i.e., it seems that the lines are broader in the upper 
parts of the loops than near their feet. This contrasts the results 
of the single Gaussian fit. At first sight this seems contradictory; 
however, this is due to the minor component generally show- 
ing stronger blueshifts in the footpoint regions. The line width 
and its implication for the coronal heating and dynamics are dis- 
cussed in detail in Sect.P 



3.3. Free double Gaussian fits 

Double Gaussian fits have been performed to quantify the con- 
tribution of the wing excess. First a free fit is described before 
a fit for a restricted parameter range is discussed in the follow- 
ing section. Actually, for the genetic algoritm used here for the 
fit (cf. Sect. [2]) , a range of parameters has to be specified. For 
the free fit, the parameter range was chosen broad enough not to 
constrain the fit in any way. 

Without constraints, the line profile is typically decomposed 
into a major narrow component accounting for the line core and a 
minor component accounting mostly for the wing excess, since it 
is much broader than the core component. Examples for the free 
double Gaussian fits can be found in Fig.[3](top row) for average 
spectra in about 10"xl0" regions and in Fig. [4] (top row) for 
individual spectra recored in one spatial pixel («l"x 1"). 

In general, the free double Gaussian fits are much better than 
the single Gaussian fits discussed above, with the x? on average 
a factor a three smaller (Fig. [3 left panel). The minor component 
typically contributes about 10% to 20% to the total line emission 
(right panel). Inspection of the individual fits to the spectra, of 
which here only four are shown (Fig.lU top row), shows that the 
free double Gaussian fit indeed provides a good model to fit the 
line profile (cf. Sect. 13.51) . 

The spatial maps for line shift and width for both compo- 
nents are shown in Fig.|2]in the middle and right columns. Given 
that the minor component contributes less than 1/5 to the line, 
one would expect the line shifts of the core component for the 
double Gaussian fit and the shift of the single Gaussian fit to 



3.4. Forcing a high-velocity component 

The double Gaussian fits reveal a high-blueshifted component in 
restricted areas at the periphery of the active region, while in the 
central part covering most of the area the second minor compo- 
nent shows moderate blueshifts combined with large line widths. 
Considering a spectrum as shown in column b of Fig. (3) one 
could still argue that a double Gaussian fit with a high-blueshift 
component would do as good a job. Thus a double Gaussian fit 
was performed with two constraints: (1) The minor component 
has to be shifted by more than 60km/s to the blue and (2) its 
Gaussian width should be smaller than 65 km/s (i.e. not larger 
than the maximum width of the core). 

The resulting Doppler shifts of the minor high- velocity com- 
ponent are in the range of 100 km/s to 140 km/s: if the minor 
component is forbidden to also account for the excess in the 
red wing, then it moves far out to the blue. This is depicted in 
the left panel of Fig. [6] showing the map of shifts of the minor 
high-velocity component. Very high blueshifts on the order of 
100 km/s are found everywhere. That the shifts are always be- 
low 160 km/s confirms that the line profile of Fe XV investigated 
here is not blended by the Fe xvn line (cf. Sect.0. As Fe xvn 
would be the strongest potential blend, it can be concluded that 
the present study is not affected by line blends. 

When looking at individual spectra, it quickly becomes clear 
that the fit constrained to have a high- velocity component is not 

3 Some of the spectra i n the se areas might be affected by line bends, 
which is discussed in Sect. l4.1l and in more detail in appendix |A| 
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Fig. 6. Constrained double Gaussian fits and 
their significance. The left panel shows the 
Doppler map for the minor (high-velocity) 
component of the Gaussian fits constrained to 
have a high- velocity component. 
The right panel shows contour lines of the ra- 
tio R x of the Xr f° r me constrained double 
Gaussian fit divided by the Xr f° r me f ree dou- 
ble Gaussian fit. Same field of view as in Figs.[T] 
andEl See Sect.EH 



doing very well in general. It simply fails to properly account 
for the red wing emission (see panels f31, 7I, SI in Fig.©. Only 
in those restricted regions where the free double Gaussian fit 
shows a strongly blueshifted component does the constrained fit 
do a good job (Fig.ffl column a). This is underlined by the x 2 
values given in the panels of Fig.lU as well as by the histogram 
for the Xr - m general (outside the high-blueshift regions), the x 2 
for the constrained fit is about a factor of 2 larger than for the 
free double Gaussian fit. This is even true when looking at the 
footpoint regions of the hot loops, e.g. in region (b) outlined in 
Figs.[Uand[2 

To emphasize that the major part of the active region is bet- 
ter fit by a free double Gaussian, contour lines for the ratio R x 
of the Xr f° r m e constrained fit divided by the one for the free 
double Gaussian fit are plotted in Fig. [6] Inside the red contour 
line for R x =l.2, the free double Gaussian fit is certainly more 
significant than the fit with the constrained high-blueshift com- 
ponent, and this region covers the whole part of the active region 
with bright coronal emission (cf. lower right panel of Fig.Q]). 
Outside the R x =l.O contour, where the constrained fit is better 
than the free fit, only the high-velocity regions are found (thick 
black contours). These are enclosed by a type of pocket of the 
R x =l.O contour. (R x is also below 1 in areas with very low sig- 
nal in the coronal line, where the fit is not reliable for both the 
constrained and the free Gaussian fits.) 

From this one can conclude that a high-blueshift compo- 
nent is indeed only present in the two restricted areas outlined 
in Figs.Q]and[2]by the contour lines, which are located at the 
periphery of the active region. In the major part of the active 
region, the minor component of the double Gaussian is broad 
(more than twice the line width of the core) and shows moderate 
blueshifts ranging from zero to some 50 km/s (cf. Fig. [8]). 

3.5. What is the best fit to the line profile? 

In the preceding sections, the reduced Xr was use d to 
judge the significance of th e fit. This is defined as (e.g. 
Bevi ngton & Robinsonl Il992h 



N 



1 df — m(vi) ' 



N-f 



(3) 



Here di denotes the data (i.e. spectral radiance) taken at spec- 
tral position i and m(vi) is the model value of the fit at the 
Doppler shift V{ of that position, where m is a (single or dou- 
ble) Gaussian as defined in (Q} plus a continuum. The respective 



measurement errors are given by oi . Data at 7V=24 spectral posi- 
tions have been used and the fit has f=7 (4) degrees of freedom 
for a double (single) Gaussian fit. Assuming that the individual 
measurements (i.e. the spectra of the spatial scan) are indepen- 
dent (and the errors are normally distributed), in a large enough 
data set the values of x 2 s hould be distributed with a pea k around 
(for basics see e.g. lBevington & Robinsonl fl 992). 

The Xr distribution for the best fit, i.e. the free double 
Gaussian, shows a peak near and roughly f ollows 

the theoretical distribution (iBevington & Robinsonl Il992l Sect. 
11.1) if scaled by a factor of 2 (cf. left panel of Fig.0. It 
is unclear why exactly the peak of the Xr distribution is only 
0.5, i.e. less than 1. One reason could be that the errors <Ji 
of the data are overestimated. The errors returned by the EIS 
data calibration software are basically Poisson errors, but they 
also account for o ther effects such as CCD read-out noise (see 
I Young et aill2009l) . If these errors were overestimated by some 
40% (p^l—y/l, this would account for the factor of 2 in the x 2 
distribution (cf. Eq.0. A detailed analysis of the measurement 
errors of EIS would be needed to fully investigate this, which 
is beyond the scope of this paper. The assumption of a normal 
distribution of the errors for the x 2 test is a potential problem, 
too. It remains open as to what extent this will influence the dis- 
tribution of x 2 and in which direction. Possible problems in the 
determination of the measurement error and for the underlying 
assumptions for the x 2 test show that it is no surprise that the 
peak of the x 2 distribution is not at 1 . 

Despite these problems, one can compare the x 2 distribu- 
tions for the different fits. From Fig. [5] (left panel) it is clear, 
that the free double Gaussian performs best, while the single 
Gaussian fit clearly fails to be a good representation of the data 
since it cannot account for the wing excess. The double Gaussian 
fits constrained to have a high-blueshift component perform bet- 
ter than the single Gaussian fits, but because the high-velocity 
component is only present in (two) restricted areas, the con- 
strained Gaussian fit falls short compared to the fit with two free 
Gaussians. 



3.6. Errors for best-fit parameters 

When using the genetic algorithm-based optimization one sould 
have to investigate the curvature of the x 2 surface in the multidi- 
mensional parameter space to determine the errors of the line-fit 
parameters. As this is a complex problem of its own, a local 
optimization method is used with the results from the genetic al- 
gorithm as the initial guess to estimate the uncertainties of the fit 
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Table 2. Contribution of possible blends to Fe XV at 284 A. 



rest 


relative 


modeled contribution [%] a 


wavelength 


shift 6 




active 


quiet 


coronal 


line [A] 


[km/s] 


flare 


region 


Sun 


hole 


Fexvil 284.010 c 


-158 


2.8 


0.1 


<0.1 


0.0 


Mix 284.015 d 


-152 


0.2 


0.6 


2.1 


83.7 


Fexv 284.160 6 





97.0 


99.3 


97.8 


16.3 



a Contribution modeled using the Chianti atomic data base. 

b Doppler shift according to difference in rest wavelength. 

c From NIST atomic spectra data base (Ralch enko et all 2008). 

d From NIST (as c ); based on laboratory measurements by 

Val ero & Goorvitchl d!972h . who estimate uncertainty to be ±0.04 A 

(±40 km/s). 

e Based on solar observations as provided bv lBrown et all d2008h . un- 
certainty is given as ±0.004 A (±4 km/s). Within errors consistent with 
NIST data base (observed and theoretical calculated values). 

(cf. Sect. [3 . This should provide a good approximation for the 
errors of the fit. 

In the case of the free double Gaussian fit, the mean errors for 
the line shift are ±5 km/s for the core and ±7 km/s for the minor 
component. The respective values for the errors of the Gaussian 
width are ±7 km/s and ±16 km/s. It is clear that the errors are 
greater for the minor component, since the shift and the width 
of a weak (mostly broad) component are more difficult to deter- 
mine. Please note that these are just the errors for the line profile 
fit so they do not include other error sources. 

4. Discussion 

4. 1. Possible effects of line blends 

Acco rding to the Chianti atomic data base (v6.0; iDere et all 
I2009I) . there are two potential blends in the range of 300 km/s 
to the red or blue of the Fexv line investigated here: an 
Fexvil line fo r med in hot plasmas such as found in flares 
(Feldman etall Il985t IPeT Zanna & Ishikawa, 2009) and an 
Al IX line formed at som ewhat lower temperatures than Fe XV 
(e.g. jBrown et all [2008b . To estimate the potential influence of 
these blends the relative contribution of the three lines was cal- 
culated using Chianti for different situations (using different dif- 
ferential emission measure curves). The results are summarized 
in Table[2] The Fe XVII line can be ruled out since it would play 
a role only in a flare, which is not considered here. Under ac- 
tive region conditions, the Al IX line should contribute less than 
a percent, and even in the quiet Sun, it sho uld be a factor of 
50 weaker than Fexv. As discussed in Sect. l3.3l about 10% to 
20% of the line intensity is excess emission in the line wings 
(cf. Fig. right panel for free double Gaussian), and the excess 
emission shows shifts well below the approx. 150 km/s shift of 
the potential blend (cf. Fig.0 Fig.[7J and Table©. 

Thus it can be concluded that in the parts of the active region 
which are bright in Fe XV the blending by Al IX is negligible. 
Any deviations from a single Gaussian line profile which are 
larger than about 5% cannot be attributed to blends, but have to 
result from flows along the line of sight. However, in the areas 
with low emission in Fe XV, which might be relatively cool, the 
blending by Al IX could be present. Going to the extreme case 
of a coronal hole with maximum temperatures below 1 MK, the 
Al IX emission would even dominate Fe XV (cf. Table©. In the 
regions where the minor component of the free double Gaussian 
fit shows high blueshifts at low Fe XV emission, in about 1/3 of 
the area the fit might pick up the blend. This is discussed in detail 



in appendix [Aj The is sue o f a potential Al IX blend affects only 
the discussion in Sect. 14.41 and then only partly, as the blend is 
only found in a fraction of the region discussed there. 

4.2. Instrumental effects causing line asymmetries? 

If an asymmetry of the instrumental profile caused the excess in 
the line wings, this effect should be roughly the same across the 
field of view. However, the broad minor components have quite 
different strengths, shifts, and widths. For example, in the aver- 
age spectra shown in panel bl, cl, and dl in Fig.O the relative 
shift between the core and the minor component spans a range 
from 7 km/s to 35 km/s. This is considerable compared to the 
width of the line core and on the same order as the (exponential) 
width of the instrumental profile (assuming it were Gaussian). 
This wide range of relative shifts between minor and core com- 
ponents is a clear indication that the effect is non-instrumental, 
especially as the different relative shifts are associated with solar 
structures, i.e., loop footpoints and loop tops (cf. panels bl and 
cl in Fig.©. 

4.3. Asymmetric velocity distribution function of the ions? 

An asymmetric velocity distribution function of the emitting 
ions could be the reason for the line asymmetries, too. Usually 
it is assumed that the high densities in an active region ensure 
strong enough collisional coupling to have Maxwellian distribu- 
tion functions. How ever, deviations from this seem to be possi- 
ble, as shown e.g. by lBourouaine et al"] (|2008a) for non-resonant 
ion heating by Alfven waves in the low corona. The extended 
wings of the ion distribution in velocity space would directly 
translate into the observable emission line profile. Were this the 
cause of the asymmetric line profile, the interpretation would be 
radically different from the interpretation discussed here. Instead 
of spatially separated source regions for the line core and the 
minor component (cf. Fig.fTOl), the emission could originate in 
a monolithic structure. The asymmetry would not arise from 
structures in real space, but from peculiarities in velocity space. 
Future investigations will have to show whether this process is 
relevant in active regions, too, and if it indeed is compatible with 
the observations presented here. 

4.4. High-velocity upflows in the active region periphery 

The minor components of the double Gaussian fits fall into two 
categories. (1) In the high-emission part of the field of view, i.e., 
in the core of the active region, the minor component is broad 
and shows moderate shifts up to some 50 km/s, which is dis- 
cussed in Sects. l4.5l and l4.6l (2) Distinct from this are two small 
regions at the periphery of the active region with narrower mi- 
nor components at high blueshifts of about 40 km/s to 60 km/s, 
which is discussed in this section. There is also a fraction of area 
(about 1/3) in these regions that shows shifts of some 100 km/s, 
but it is not clear that this is truly a component of Fe X V or a 
blend (see discussion in appendix [A] and histogram in Fig. IA.2l) . 

Both high- velocity regions are located in areas with low 
coronal emission at 2.5 MK (logT=6.4) and are outlined by con- 
tour lines near (a) and (a) in Figs.Q]and[2] It is important to note 
that the line core emission, accounting for about 80% to 90% of 
the total line emission, shows very low blueshifts, or even red- 
shifts in these high- velocity areas (cf. Fig. [2 middle and right 
panels in top row). In both cases one might consider the line 
core component being emitted in the major part of the structure, 
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where plasma is sitting in a more or less static fashion. As seen 
in Fig.Q] both areas are outside the photospheric magnetic field 
concentrations, and they do not seem to be related to the emis- 
sion from the chromosphere and transition region seen in He II. 
However, they are quite different in terms of emission from the 
cooler parts of the corona. 

4.4.1 . Solar wind outflow from filamentary funnels? 

After comparison with the FeVlll line formed at 0.5 MK 
(logT=5.7) in Fig.Q] the high- velocity region near (a) seems to 
be above a fan-like cool structure. This might be interpreted as 
a coronal funnel, the cool emission indicating the location of its 
footpoint, but further analysis of magnetic data would be needed 
for a conclusion on this. Even though not emerging from a full 
sunspot, this structur e might also be related to a sunspot plume 
(Maltbvetal., 1999) because the emission from Fe x shows a 
redshift associated with it (cf. Fig. IA.lb . 

Since the line core emission of Fe XV from the major part 
of the structure is not shifted (or only little), the high- velocity 
component could then be considered as "fingers" reaching into 
the main structure. These fingers could either supply mass to 
the structure or they could be the base of strands through which 
material is channeled into the upper corona or even the wind. 
Especially the high- velocity region near (a) might be a good can- 
didate to supply mass to the wind at high velocities already close 
to the surface. The temperature in solar wind streams should be 
rather low ( <1 MK) as show n for fast streams from polar coro- 
nal holes by IWilhelm et al.l (Il998l) . because most energy goes 
into acceleration instead of heating. Thus the very low Fe XV in- 
tensity in this high- velocity region near (a) as seen in Fig. Q] or 
IB. H is in line with the interpretation that this represents the origin 
of a solar wind stream. 

In contrast to previous scenarios, here the wind outflow 
would come from small strands that are embedded in an almost 
static coronal funnel. This would cha nge the pictur e of a more 
or less homogeneous coronal funnel ( Gabriel, 1976) at the base 
of the wind to a scenario where the funnel would be a filamen- 
tary structure with fast outflows and more static structures next 
to each other. In a way one could look at this as a smaller version 
of plume and inter-plume regions in (polar) coronal holes where 
the origin of the fast wind is found. 

It should be noted that the findings presented h ere are quite 
differ ent from the studies of iDoschek et al.l (120081) or iHe et aTl 
( 2010) who also find evidence of a wind outflow in the periph- 
ery of an active region. However, they only investigate single 
Gaussian fits. The regions with high- velocity components in the 
double Gaussian fits such as are found here coincide with the lo- 
cations of strong blueshifts of the single Gaussian fits (cf. Fig. [2]) 
— the single Gaussians show the strong blueshifts because the 
optimization procedure does its best to fit the line excess in 
the blue wing and thus leads to a strong blueshift of the single 
Gaussia n. Thus one can speculate that the lDoschek et al.l ( 2008) 
and the lHeetal.1 (1201(1) data actually show the same effect of 
a high-velocity component, which will have to be checked by 
re-examining their data. 

4.4.2. Nanoflare heated loops? 

The other high- velocity region located near (a) in Fig.Q] is lo- 
cated above strong emission in Fe VIII, as well as in Fe X formed 
at 1.2 MK (logT=6.1). One might speculate wheather this high- 
velocity region is above a footpoint of a faint hot loop following 



the crosses (a) - (5). It is barely visible in Fe XV when plotted on 
a linear scale ( Fig.Q ]), but clearer when displayed on a logarith- 
mic scale (Fig. lBJl) . The strong emission at lower temperatures 
might hint at an increased level of heating at the base of the loop, 
e.g. in the form of field line braiding tha t would lead to high en- 
ergy deposition low in the atmosphere dGudiksen & Nordlundi 
120021) . 

For an increased activity near the loop foot, one can expect 
that the braiding and twisting of the field lines will result in fre- 
quent nanoflare events in the coronal part of the loop. The ob- 
served high-blueshift component near the loop footpoint would 
be in l ine with the model prediction by lPatsourakos & K limchuk 
(2006) based on a model for a nanoflare heated loop. In this 
model a high- velocity upflow at the footpoint (i.e. small chromo- 
spheric evaporation) is driven by the nanoflare-induced heat in- 
put in the coronal part of the loop. While the relative strength of 
the spectral component due to the high- velocity upflow in their 
model is similar to what is found in the observations presented 
here, the Doppler shift in the model is about twice as high. As for 
the filamentary funnels, here one would expect that the upflow 
would happen in a (small) part of the individual strands form- 
ing the hot coronal loop, which again would be in line with the 
iPatsourakos & Klimchu k (2006) model. 

Thus one could conclude that these observations sup- 
port flie_jianojlare_Jie^ scenario as proposed 
by IPatsourakos & Klimchukl (120061) (assuming that the high- 
velocity component indeed is not due to a blend). However, 
there are also problems. These high-blueshift regions are seen 
in a small part of the field of view and only at the periph- 
ery o f the active region, which would challenge the uniqueness 
of the IPatsourakos & Klimchu k ( 2006) model. Furthermore, the 
loop with the high-upflow leg is quite dim, while the heart 
of the active region with much stronger emission from hot 
plasma does not show the pronounced high-velocity compo- 
nent s. Thus it might be questionab le if the process described 
by IPatsourakos & Klimchukl (120061) were the dominant heating 
mechanism. 

However, it should be emphasized that, in principle, the high- 
velocity components as predicted by Patsourakos & Klimchuk 
( 2006) could still be hidden below the intense radiation from the 
main part of the active region. More precisely, it could be that 
the high- velocity component is not visible because other effects 
causing the excess in the line wings are stronger as discussed 
in the next section. Furthermore, their model is a pure (M)HD 
description of flows along the loop. In a more realistic model, it 
might be that they will find a stronger broadening of the emission 
from the upflows. For example, some MHD turbulence induced 
by the evaporation process could lead to a line broadening and 
might decelerate the upflow, which could be the reason for the 
observed properties of the broad minor component in the bight 
active region loops. 

Thus the present observations cannot be used to rule out 
their attractive proposal. Future work will have to show whether 
the upflows induced by the nanoflares could be identical to the 
source region of the broad minor components in the active re- 
gions loops. 

4.5. Heating and mass supply in the active region 

For the remainder of the paper, the discussion turns to the major 
part of the active region bright in Fe XV, where most of the coro- 
nal emission originates. T his r egion is not susceptible to effects 
of line blending (cf. Sect. l4.1l and appendix lAl). Here the spec- 
tra are composed of a narrow line core and a minor component. 
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Fig. 7. Distribution of line shift and line width 
for free double Gaussian fits to Fe XV (284 A). 
The distributions are shown for the line core 
component as well as for the minor component 
labeled "wing". The minor component has a 
contribution at high blueshifts, as can be seen 
in the inlet in the left panel, which is most prob- 
ably due to a blend (cf. appendix |A]). The small 
distinct contribution to the line widths at very 
small widths (40km/s) is an artifact of bad fits. 
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Typically, the latter contributes 10% to 20% to the line intensity, 
is a factor of two borader than the core, and shows line shifts 
from to some 50 km/s to the blue, while the core component 
shows only little shifts lower than ±10 km/s (cf. Figs. [5] and [7]). 

4.5.1 . Waves and velocity anisotropies 

Both the core and the minor component show a trend toward 
weaker line broadening near the footpoints of the coronal struc- 
tures (see Sect. l3.3b . This could either the stronger broadening 
effect higher up in the coronal loops at lower densities, or the 
broadening could be due to an anisotropic velocity distribution 
of the ions. 

Assuming a constant (or only slowly varying) energy flux 
density, a compressible disturbance (e.g. magnetoacoustic wave) 
moving upwards along a loop would steepen in amplitude as 
the density decreases. This is well known in the chromosphere, 
where acou stic waves steepen into shocks and form the Ca II 
grains (e.g. ICarlsson & SteinL [l997). A similar effect could in 
principle be expected for the corona. In the framework of coro- 
nal seismology, models and observations (mostly from imaging 
instruments) are used to investigate these periodic disturbances 
in the corona that h ave recently also been found in complex 3D 
configurations (e.g. [Of man, 200^). 

Longitu d inal w av es such as discussed recently by, e.g, 
iDe Moortell (l2QQ9h or iDe Moortel & Bradshawl (I2008I) do not 
seem to be important here, because they should produce a 
stronger line broadening near the loop footpoints, where the an- 
gle between the magnetic field and the observer is relatively 
small. Higher up, where the line of sight is perpendicular to the 
loop magnetic field, the broadening effect of a longitudinal wave 
should be very small. Thus the predicted change of line broad- 
ening would be opposite to what is observed here. 

In conseque nce, transverse magnetoa coustic oscillations 
(e.g., review of iRuderman & ErdeTyil l2009l) would predict just 
the observed trend from loop foot to apex. However, in the mod- 
els these waves are damped to heat the coronal plasma. Thus 
their amplitude should decrease with height in the loop, which 
would result in a rather small line broadening. Further modeling 
efforts would be needed to investigate what the models predict 
in terms of line broadening along a loop-like structure, which 
wou ld only imply small extensions of the modelin g effor ts of 
e.g.Uerradas et al. ( 2008) orl Van D oorsselaer e et a 1. (2009). 

Alfven waves in the corona have been directly detected 
through coronagraph ob s ervati ons of the magnetic field only 
recently dTomczyk et all 120071) . This was the first proof that 
Alfveninc disturbances propagate upward and downward along 
loop structures. As these are transverse waves, they would nicely 



explain why the broadening is the smallest at the footpoints 
where the line of sight is (more or less) parallel to the mag- 
netic field and thus perpe ndicular to the veloci ty disturbances. 
For transition region lines, Mcint osh et ID (2007) present a new 
interpretation for the observed line widths that relies on such 
Alfvenic motions superposed by longitudinal motions along the 
line of sight. Alfvenic motions have been foun d for the long- 
known type I spicules (De Pon tieu et all 12007 ah. as well as for 
the mo re recently identified type II spicules (De Pontieu et ID, 
200 7b|). However, th e coronal Alfven waves as identified by 
iTomczyk et al.l (120071) have periods of some 5 minutes and thus 
would not cause a line broadening of the coronal spectra investi- 
gated here, which have 30 s exposure time. (For a multi- stranded 
loop such low-frequency waves might have an effect even for 
shorter exposure times if many strands are in one resolution el- 
ement of the detector.) It remains to be seen whether Alfvenic 
disturbances with high frequencies and high enough power exist 
in the corona to explain the observed line broadening. 

Kinetic models for heating the corona and accelerating the 
wind by high-frequency Alfven waves through ion-cyclotron 
resonances predict high temperatures of the ions, in particular, 
the temperature perpendicular to the magnetic field would be 
considerably higher than those parallel to the magnetic field, as 
the velocity distribution functions will be anisotropic. This has 
been discussed in de tai l for the open co rona extensively e.g. by 
iMarsch & Tul (1 19971) or lCranmed (120001) . and evidence for the di- 
rect observation of the broadening throu gh this process close t o 
the Sun above the limb was presented by lPeter & Vocksl (120031) . 
More recently, models for resonant and non-resonant wave- 
particle interaction showed that also in closed large-scale mag- 
netic structures str ong temperature anisotropi es can be present 
for the heavy ions dBourouaine et al.l.[2008allbl) . The anisotropic 
velocity distribution in this kind of models with larger (average) 
velocities perpendicular to the magnetic field would be consis- 
tent with the line broadening being larger towards the loop top. 

As for the Alfven waves discussed above, further modeling 
would be needed to investigate if line broadening predicted by 
this anisotropy would agree also quantitatively with the observa- 
tions, and if the kinetic models are indeed applicable in the high- 
density environment of an active region as investigated here. 

4.5.2. Mass supply to coronal structures 

A central question, which has not been addressed yet, con- 
cerns the nature of the line core and the minor component. 
Recen tly, IDe Pontieu et al.1 (l2009h and iMcIntosh & De Pontieul 
( 2009) have looked at the asymmetry of line profiles by compar- 
ing the emission from the blue and the red wings at the same dis- 
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tance from the line centroid (integrating over some wavelength 
range). Through this they find and confirm the asymmetry as 
found by Hara et al. (2008) in other data sets and adopt the inter- 
pretation of lHara etaD (120081) that the line asymmetry is caused 
by a high- velocity upflow in the blue wing. They associate these 
strong upflows with type II spicules. Their interpretation would 
be consistent with the results shown here, with some important 
alterations and extensions. These are mainly based on the more 
detailed data analysis, which allows also quantitative measures 
for line shifts and widths. 

The upflow speeds of 50km/s to 100 km/s at the loop foot- 
points that are derived by iDe Pontieu et al ] (120091) are signifi- 
cantly higher than the blueshifts presented in this paper for the 
footpoint regions. This difference can be reconciled by inspec- 
tion of column (b) of Fig. [3] showing an average spectrum of a 
lO^xlO" footpoint region. In that example the proper double 
Gaussian fit returns a minor component with a blueshift of about 
50 km/s. In contrast first fitting a single Gaussian and then inves- 
tigating the residual (panel b3) would indicate an upflow with 
more than 100 km/s (with a very small line width). This dif- 
ference is simply because the single Gaussian does its best to 
also account for the wing excess and thus overestimates the con- 
tribution from the line core. Already assuming a high-velocity 
component will give high- velocity narrow minor components, 
which are less significant than the double Gaussians with a free 
(broad) minor component (cf. Sect. l3.4b . Comparing the typical 
Doppler shifts of the minor component near the footpoints of 
some 30 km/s to 50 km/s one can conclude that lDe Pontieu et all 
(120091) might overestimate the upflow speed by a factor of 2. 

More importantly, the present study can give information on 
the line width of the minor component, and thus provides a way 
to investigate the physical mechanism driving the mass upflow 
as discussed above in Sect. l4.5.ll 

Future studies comparing and combining different meth- 
ods to characterize the line a symmetries have to be con- 
ducted. The method used b y IDe Pontieu et al.l d2QQ9h and 
Mcintos h & De Pontie u ( 2009) is quick and very robust, espe- 
cially if the data become noisy. The study presented here re- 
quires clean spectral profiles to give reliable results but returns 
more quantitative information. 

4.5.3. Statistical connection between line shift and 
broadening 

The relation of line shifts and widths contains valuable informa- 
tion on the heating and dynamics of the plasma and might reveal 
details of the governing processes. For example a clear relation 
of increased line width with high Doppler shift could indicate a 
heating process driving a flow. 

When plotting a two-dimensional histogram of line shifts 
and widths for the line core of the free double Gaussian fit, no 
clear relation can be found (Fig.[!i)0 One might argue that there 
is a weak relation, and a closer inspection shows that the me- 
dian width of regions with positive shifts (upflows) is ~49 km/s, 
which is a little larger than the median width of ~47 km/s in 
redshifted (downflow) areas. This weak relation of shift to width 
indicates that there is no clear connection between heating of 



4 It should be noted that for a single Gaussian fit a clear relation is 
found (with higher blueshifts at larger widths; not shown in this paper). 
This is only because the line wing excess is not properly accounted for, 
and thus is an artifact of the fit. It is related to the increased line widths 
and shi fts s een for the single Gaussians in the footpoint areas (cf. end 
ofSect.HU. 
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Fig. 8. Two-dimensional histograms of line width vs. shift for 
free double Gaussian fit to Fexv (284 A). The top panel (a) 
shows the line core component, the bottom panel (b) the minor 
or wing component. The histograms are plotted on a logarith- 
mic scale, the contours showing levels of 0.1 (dashed) and 0.03 
(solid) of the maximum value. For the minor component a con- 
tribution at high shifts with small widths most probably is due to 
a blend (cf. appendix [A}. The peak with zero shift at very small 
width (lo wer rig ht part o f panel b) is an artifact due to bad fits. 
See Sects. l4331 andl431 



the plasma and the associated flows in the regions where the line 
core of the emission profile is formed. 

The situation is different for the minor component in the 
major part of the active region. Besides a large part of spectra 
with very low shifts (lower ±10 km/s), there is a tendency for 
the minor component to have a higher blueshift for increasing 
line widths (see Fig.[8]3, solid contour line). This indicates that 
a stronger heating in the source region of the minor component 
leads to increased line broadening and drives an upward flow 
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into the corona. This interpretation suggests that the actual mass 
supply to the corona is located in the source region of the minor 
compone nt. This interpreta tion would be consistent with the pro- 
posal by iDe Pontieu et al.1 d2009h and iMcIntosh & De Pontieul 
( 2009|) if the type II spicules they propose as the mass suppli- 
ers were associated with the broad minor spectral component of 
the coro nal lines. As these type II spicules host strong Alfven 
waves (IDe Pontieu et al.[|2Q0 7b). this seems reasonable. 

4.6. Self-similarity from the transition region to the corona 

An active region as investigated in this paper and a quiet Sun net- 
work patch at transition region temperatures are quite different 
in terms of spatial structure, size, and energetics. Nonetheless, 
numerous similarities of the spectral profile properties (in a sta- 
tistical sense) can be identified. This hints at common processes 
on very different scales. 

For transition region lines iPeterl (I2Q0QL 1200 ll) analyzed the 
spectral profiles in a similar way to the present paper, i.e., dou- 
ble Gaussian fits were performed. In network regions, these also 
typically show a narrow line core and a broad minor component. 
The minor component contributes some 10% to 20% to the total 
intensity and is about a factor of 2 to 3 broader than the core 
component, consistent with the coronal line investigated in this 
paper. 

For the Civ line at 1548 A formed at around 0.1 MK, the 
relation of widt h to shift for both spectral components has been 
investigated by IPeterl (l2Q0Qh . The absolute values of shift and 
width are lower, but their relation is remarkably similar to the 
coronal Fe XV line: while no relation is found for the line core, 
the line shift is closely related to the line width, with low shifts 
for small widths. 

Finally, the line width for the minor component found for 
the coronal line in this st udy fa l ls into the trend found with 
the transition region lines. Peter] (120011) finds a monotonic in- 
crease in the non-thermal line widths for the minor component 
of the transition region lines (his Fig. 7) following a power law 
T 1 / 4 (with the line formation temperature T). Extrapolating to 
logT=6.4 where Fexv is formed, this gives about a lOOkm/s 
non-thermal speed. At the footpoint regions of the active re- 
gion investigated here with blueshifts of the minor component 
of 40km/s to 50km/s, the exponential width of the minor com- 
ponent of Fexv is found to be about HOkm/s to 120km/s (cf. 
Fig -Eh)- Using the instrumental and thermal widths as listed 
in Table[TJ the non-thermal speed of Fe XV is in the range of 
lOOkm/s to HOkm/s, just matching the extrapolation from the 
transition region. 

There is one major difference, though. While the medium 
value for the line shift of the minor components in the transi- 
tion region is low (and more or less consistent with zero shift, 
considering the error bars), the minor components of the coronal 
line show a clear blueshift of up to 50 km/s (which is about 20% 
of the sound speed at logT=6.4). This could indicate that the ac- 
tual mass feeding to the coron a is located at te mperatures above 
several 10 5 K, as suggested by lTu et all (120051) . 

The above similarities of the statistical properties and the 
match of the (coronal) line width to the extrapolation from the 
transition region suggests that there is a common process gov- 
erning the physics leading to the minor components of the line 
profiles all the way from temperatures of 10 4 K to several 10 6 K. 
In a self-similar fashion, network patches in the transition region 
and large active regions seem to operate in a unique way. This is 
remarkable because the magnetic flux in a single network patch 



is many orders of magnitude smaller than in an active region. 
Future work on more coronal lines formed at various tempera- 
tures and on different active regions will be needed to solidify 
this picture. 

4. 7. Flows in loops matching 3D MHD model 

The observations presented here should also be discussed briefly 
in the context of 3D MHD models of the coronal structure and 
dynamics. This discussion is to show the potential of combin- 
ing these observations and models. Further investigations will 
be needed to fully explore this. 

The emissi on from such a 3D MHD coronal model has been 
synthesized by lPeter et al.l d2004. 2006). This is based on heating 
the corona through the dissipation of currents induced by braid- 
ing the (coro nal) magnetic field through photosph eric connec- 
tive motions dGudiksen & NordlundL I2002L l2005allbh . As these 
3D coronal models cannot (yet) resolve the small structures that 
might host the source region of the minor components (e.g. 
type II spicules), the discussion concentrates on a brief com- 
parison between the shifts seen in the line core component to 
those derived from the MHD model. In Fig.[9]the line intensity 
and shift synthesized for the Mg X line are displayed with the 
computational box viewed from above. The data are taken from 
iPeter et aLl (|2006), who also provides movies of the temporal 
evolution. 

In the observation, there are two distinct flow patterns to be 
found within the same active region. The hot loop system seen 
in the Fe XV line in Fig. [T] with the apex close to (c) shows a nice 
siphon flow (see middle top panel in Fig. [3. This flow is proba- 
bly driven by asymmetric heating at the loop feet, and a further 
investigation would be needed to see whether this is really sup- 
ported by the line width observations. A hot faint loop can be 
seen in Fe XV with its apex close to 7 (more clearly seen with 
logarithmic scaling in Fig. IB.lb . This shows downflows on both 
sides of the loop, indicating that the loop drains, probably be- 
cause it is cooling down (cf. different loop clases identified by 
lUgarte-Urra et al.L l2009h . 

Such loop structures can be modeled in the framework of ID 
loop m odels with a prescrib ed heating function. For example, 
recently IWarren et all (120101) has presented a class of ID multi- 
thread models in which they successfully reproduce observations 
by adjusting the energy release time (and by this prescribing the 
heat input). In contrast to this, the temporal and spatial distribu- 
tions of the energy input, as well as the resulting spatial struc- 
tures, are calculated self-consistently in the 3D MHD model. 
(For example, at the cost of lower spatial resolution, which is 
also affecting the thermal conduction; the ID and 3D models are 
complementary.) 

In the 3D MHD models, the same flow patterns can be found 
as are seen in t he active re g ion s tudied here. In the synthetic 
spectral maps of IPeter et alJ (120061) . one can find siphon flows, 
as well as draining loops. Examples roughly matching the ob- 
served cases are highlighted in Fig. [9] and labeled (c) and (7) 
accordingly. While these structures are present in the model at 
the time of the snapshot shown in Fig.O the situation is quite 
different only 10 m inutes later. This would be consistent with 
the observations by lUgarte-Urra et aLl (l20Q9h . 

This comparison of a 3D MHD model snapshot to observa- 
tions shows that the 3D MHD models can in principle repro- 
duce the flow patterns seen in real active regions. However, the 
magnetic configuration at the lower boundary of the 3D MHD 
model and the photospheric magnetic field for the coronal obser- 
vations discussed here are different. For proof, new simulations 
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Fig. 9. Comparison to re sults from 3D MHD coronal model. Displayed are the synthesized coronal emission and Doppler shifts 
(iPeter et all 120041 120061) for the Mg X line at 625 A. This snapshot of the time-dependent simulation shows the computational box 
being viewed from straight above. The intensity is normalized to its median value (J) and plotted on an inverse scale. The Doppler 
map is plotted on the same scale as the map of the observed line core shifts in Fig. [2] The dashed lines indicate the location of 
loops showing a siphon flow (c) or a downflow on both sides (7). These loops mi ght b e compared to the loop systems seen in the 
observations with their apex close to (c) and (7) in Figs. [Tl (or IB II) and[2j See Sect. 14.71 



will thus be needed. In these, the evolving magnetic field as ob- 
served in the photosphere will be used as a (time-dependent) 
lower boundary condition. The resulting synthetic spectra from 
the 3D model can then be compared directly to the actually 
observed coronal emission. The recently launched space-based 
Solar Dynamics Observatory (SDO) will provide an ideal plat- 
form for such investigations closely tying observations to 3D 
MHD coronal models. 

5. Conclusions 

This study emphasizes the importance of information hidden in 
details of the spectral line profiles, which show severe deviations 
from a single Gaussian in coronal lines formed in active regions. 
In this manuscript the interpretation is adopted that (two) differ- 
ent spatial structures contribute to the line profile and, through 
this, lead to a double Gaussian profile. The line fits show that 
the line profile is in general composed of a narrow Gaussian ac- 
counting for the line core emission and a minor broad Gaussian 
component for the excess emission in the line wings. The results 
for the line shift for the line core can differ significantly from a 
single Gaussian fit. In some areas even opposite results can be 
found when comparing the Doppler shifts for the line core and 
the single Gaussian (as in the region with fast upflow compo- 
nents, cf. Fig.[2]and Sect.l331). 

In two small areas at the periphery of the sunspot some evi- 
dence of a high upflow velocity component in the line profiles is 
found. This is located in regions of weak coronal emission, sup- 
posedly along low-density structures, which might be inte rpreted 
as the footpoint regions of the (slow) solar wind (Sect. l4.4.T1) . 
Evidence is found that these structures are highly filamentary 
with small fingers of high- velocity upflows feeding the wind 
sticking up into the background funnel (cf. Fig.fTOb. 

In the active region, basically all spectral profiles show a 
minor broad component that is about two times broader than 
the line core component and shows systematic blueshifts up 
to 50km/s, while the line core remains more or less unshifted 
(cf. Fig. [7]). The data indicate that the line core and the minor 
component are narrower near the loop footpoints than at the 
apex, indicating some asymmetry in the broadening mechanism. 
This could be caused by (Alfven) waves or a heating mecha- 
nism that leads to preferential perpendicular heating of the ions 
(Sect. 14311) . 



Comparing the present observations with previous work 
leads to a picture where individual strands of a loop are feed- 
ing the loop with heated material, leaving the signature of a 
broad blueshifted component in the line profile. After a while, 
these strands will no longer be heated, will slowly cool down 
and become part of the background loop that produces the nar- 
row more or less unshifted emission of the line core. Now other 
strands are heated and fill the loop, and so forth. These heated 
str ands might be identified wi th the type II spicules propose d 
bv lDe Pontieu et alJ d20Q9h and lMdntosh & De Pontieul d2QQ9h . 
This is depicted in Fig.[l0l In contrast to previous studies, this 
analysis also provides information on the heating and dynam- 
ics of the EUV counterparts of the type II spicules, which might 
be the source region of the minor components, because in the 
present study information on line widths and shifts of the minor 
component is extracted (Sect. l4.5.2t . 

Depending on the magnetic structure and the distribution 
of the heat input between the footpoints, the background loops 
might host siphon flows or they might be draining in response to 
a lower heating rate (cf. Fig.[l0|). Both cases can be found in the 
present dat a set and have also been found in 3D MHD simula- 
tions (Sect.l4~71). 

No cl ear evidence could be found to support the nano-flare 
model of iPatsourakos & Klimchukl (120061) . This model predicts 
that there is a narrow high- velocity upflow component in the 
blue wing (at some lOOkm/s, which is not present in the foot- 
point regions of the active region loops. However, further (more 
realistic) modeling efforts would be needed for any conclusive 
statement (Sect. 14.4.21) . 

The relation of line width to line shift for the line core, as 
well as for the minor component, is very similar in the coronal 
line investigated here and the transition re gion e mis sion lines 
originating in bright network patches ( S ec ts. 14331 and l4~6b . This 
shows that the basic processes driving the heating and dynamics 
on smaller scales in the network and on larger scales in the active 
region are similar in nature. (At least they have to produce the 
same spectral signatures.) This is supported by the match of the 
minor component line width of the coronal line from the active 
region to the extrapolation from the transition region lines from 
the network patch. 

In summary, these observations reveal the complex flow sys- 
tems in the active region corona leading to a mass cycle with 
small structures (probably type II spicules) feeding mass into the 
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Fig. 10. Cartoon of footpoint region of coronal structure. The left panels shows cool type II spicules (as seen in Ca II) embedded in 
the hot background plasma. These are surrounded by a thin transition region (blue) and are the basis of hot plasma flows into the 
corona (red). The strong heating and the Alfvenic motions associated with these structures lead to strongly broadened emission. This 
is then seen in the EUV spectrum as the broad minor component with strong blueshifts. The heated strands associated with the type 
II spicules feed the mass to the coronal structure. The smaller images to the right depict the coronal structures possibly emerging 
from the footpoint region (rectangles) shown to the left. The background structure near the feet can show upflows or downflows 
(depicted by the black arrows). See Sect.[5] 



corona, which could host different types of flows as depicted in 
Fig. [TO] It has to be noted, however, that this is only one possible 
interpretation. Other appealing ideas, such as an asymmetric ion 
velocity distribution function causing asymmetric line profiles 
(Sect. l4.3K will have to be followed and checked more closely 
against the new observations of coronal EUV emission line pro- 
files. 

New studies that test the ideas proposed here further will 
have to include the analysis of both emission lines with differ- 
ent mass-to-charge ratios, as well as of emission lines covering 
a range of line formation temperatures from the upper transition 
region to the hot corona. 
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Appendix A: Blending of the Fexv line 



The discussion in Sect. 14.11 made clear that, for the dataset ana- 
lyzed here, the Fe XVII line will not blend the Fe XV line of inter- 
est, as would only be the case under flare conditions. In contrast, 
the Al IX line could potentially blend Fe XV in regions where the 
plasma is cool enough. In a coronal hole, where temperatures 
stay below 1 MK the Al ix could even be the dominant emit- 
ter (cf. TableO. The present data set certainly does not include 
a coronal hole. However, in the regions where the minor com- 
ponent shows high velocities (contours close to (a) and (a) in 
Figs. [T]and Fig. [3, the temperatures will be significantly lower 
than in the main part of the active region. Nonetheless, the tem- 
peratures there can be expected to be well above 1 MK, s imply 
because there is still significant Fe XII emission (cf. Fig. IB. II) , 
which is formed around 1.5 MK under equilibrium conditions. 
According to Tabled one might estimate the strength of this 
Al IX blend to be at the level of a couple of percent. In the fol- 
lowing the discussion concentrates on these low Fe XV intensity 
high- velocity minor component regions, where the biggest effect 
of the Al IX blend can be expected. 

A major question is at what wavelength the Alix blend 
can be expected. Based on the rest wavelengths, Al ix should 
be shifted by some 150km/s to the blue relative to Fexv (cf. 
Table©. Because Al ix could show a net line shift, the Fe x line 
is used as a proxy to estimate the Al IX net Doppler shift. Both 
lines form at roughly the same temperature, so they can be ex- 
pected to have similar Doppler shifts. The Fe X line has been 
recorded co-temporally and co-spatially with the Fe XV line, and 
its line radiance is shown in Figs.[Tland lB.ll As the line profile of 
Fe X is not very clear, a single Gaussian fit was performed f or the 
line core alone. The resulting Doppler map is shown in Fig. IA.ll 
Here, as for Fe XV in the main part of this paper, the wavelength 
is calibrated by assuming an average zero shift. 

In the region of interest (outlined by the contour lines also 
in Fig. IA.llh the average shift of Fe X is some 5 km/s to the red, 
as can be seen from the histogram in Fig. IA.2l (right panel b). 
The core of Fe XV only shows a very weak shift in the region of 
interest, consistent with zero shift. Studies investigating the rel- 
ative Doppler shift between Fe X and Fe XV show that the Fe X 
line is shifted some 5 to 10 km/s to the red compar ed to Fe XV , 
depen ding on the location in the active region dTripathi et all 
2009). After combining these effects, the relative Doppler shift 
of Fe x compared to Fe XV should be to 5 km/s to the red 
(10 km/s maximum to give a margin). This should be compa- 
rable for Al IX. 

Thus if assuming a maximum relative Doppler shift of 
10 km/s combined with the 150 km/s shift of the rest wave- 
lengths of Al IX and Fe XV mentioned above, one can estimate 
that the actually observed blending Al IX line should appear at 
least some 140 km/s to the blue from the Fe XV line core. 

The histogram of the Fe XV minor component shifts in the 
high- velocity regions is shown in the left panel (a) of Fig. IA.2l 
Here a bimodal distribution is found with maxima around 
50 km/s and 115 km/s. Certainly the major part with the maxi- 
mum near 50 km/s covering about 2/3 of the area is not affected 
by the Al IX blend (even when considering the uncertainties of 
the Al IX rest wavelength of 40 km/s as given in TableO. It might 
be disputable if the part (1/3) at ^115 km/s is due to the blend 
or not. Clearly the histogram of the whole field of view shown in 
Fig-El (left panel, inlet) and by a thin blue line in Fig. IA.2l shows 
a clear peak near 130 km/s, indicating that this very minor frac- 
tion of all spectra is indeed caused by the Al IX blend. In the 
high- velocity region shown in Fig. IA.2h (thick black), one might 
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Fig. A.l. Doppler shift map of Fe x at 184 A, co-spatial and co- 
temporal with the other maps shown in this paper. The two con- 
tour lines close to (a) and (a) indicate the areas with a high- 
velocity component of the Fe XV line, just as in the other maps 
in this paper. See Fig. [T] for further details on the labeling. The 
intensit y map for this line is shown in the top right panel of Fig. Q] 
andFig.EU 



argue that a low side maximum is visible near 130 km, and thus 
the spectra with the minor component near 115 km/s might not 
be due to the blend. However, the current data do not seem to 
allow a final conclusion to be drawn on this. 

In two different active regions, iBrown et"aD (120081) find a 
component in the Fe XV spectra 114 km/s and 88 km/s to the blue 
of the Fexv rest wavelength that contribute 3% and 7% to the 
total line profile emission (their Table 2). Given the uncertain- 
ties of the rest wavelength, it could well be that these are the 
same as the part of the histogram in Fig. IA.2h close to 1 15 km/s. 
However, it seems indisputable that the high-velocity compo- 
nents found in this study for the high- velocity region centering 
on 50 km/s (being the majority) are not caused by a blend. 

It should be stressed here that, based on these estimations, 
the minor component of Fe XV in the bulk part of the active re- 
gion that is bright in Fe XV, cannot come from the Al IX line. 
The minor component has shifts that are much too low (below 
60 km/s blue) and has line widths that are too large to be due 
to a blend. The finding that the free double Gaussian with a mi- 
nor (broad) component is a better fit to the line profile th an a 
constrained Gaussian with a high- velocity component (Sect. 13.41) 
provides further evidence that in general an Al IX blend is not a 
concern in the bright part of the active region. 
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Fig.A.2. Histograms of Doppler shifts in the 
high-velocity regions outlined by the contour 
lines in the maps ; close to (a) and (a) in 
Figs.mSlOE]] The left panel (a) shows the 
histogram for the minor mostly high-velocity 
component of the Fe XV line (cf . top right panel 
in Fig. [2]), while the right panel (b) sho ws th e 
shift of the core of the Fex line (cf. Fig. IA.ll ). 
For comparison, the thin blue line in the left 
panel (a) shows the distribution of the Fexv 
minor component shift in the whole field of 
view (same as shown in Fig.[7]but scaled to 
show the contribution at high velocities). The 
vertical lines in panel (a) indicate the location 
of the maxima at high velocities near 115 km/s 
and 130 km/s of the two histograms. 



1 1 1 1 1 1 1 


R=\0(I) 
















; + 


iHe 11(256) 


log'/=4.<) 



-1.5 kG... +1.5 kG 



□ (c) 



200 

>< 

U 

o 1001 

50 



□ (b) 



□ (d) 



(a) 



+'(5) 



+(P) 



+ (7) 



r y*P\ ^ ma g net gram 



| photosphere 



-50 50 100 150 

solar X ['*] 




Fig. B.l. Same as Fig.Q] but now the intensities are plotted on a logarithmic scale. The plotting ranges R are different from Fig.Q] 



